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We present an implementation of the steady state Keldysh approach in a Green's function multiple 
scattering scheme to calculate the non-equilibrium spin density. This density is used to obtain the 
spin-transfer torque in junctions showing the magnetoresistance effect. We use our implementation 
to study the spin-transfer torque in metallic Co/Cu/Co junctions. 

PACS numbers: Valid PACS appear here 



I. INTRODUCTION 
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The discovery of the giant magnetoresistance (GMR) 
effect in metallic spin valves systems^ has led to sub- 
stantial research in the field of spintronics due to the 
possible applications including read heads in hard disks, 
storage elements in magnetic random access memory 
(MRAM), and sensors. 



An effective method for writing information into the 
elements is necessary for the application as storage el- 
ements in MRAM. In particular, one has to be able 
to change the magnetic orientation of the ferromagnetic 
leads relative to each other. One promising approach is 
the current induced switching proposed by Slonczewski 3 - 
and by Berger.— A current is driven through the junction 
and becomes spin polarized in one ferromagnetic lead. 
This polarization is conserved going through the spacer 
layer. When the corresponding angular momentum of 
the polarized current is not exactly aligned to the mag- 
netization of the second ferromagnetic lead, the electrons 
precess around the magnetic moment of the second mag- 
netic layer. In turn this precession leads to a torque 
acting on this magnetization forcing it to rotate. When 
the current is large enough, the magnetic orientation in 
the second layer can be switched. There is also a torque 
acting on the first ferromagnetic layer, but this layer is 
magnetically pinned. 



In this paper we present an ab initio calculation of the 
spin-transfer torque using a multiple scattering Green's 
function scheme. In particular, the non-equilibrium spin 
density is calculated using the steady state Keldysh ap- 
proach (see Sec. Ill Ap . This spin density is used to cal- 
culate the torque acting on the ferromagnetic layer (see 
Sec. Ill Bp . We conclude by testing our approach through 
an application to a Co/Cu/Co system, which has been 
studied by other authors^ 
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FIG. 1: Division of the junction into three regions. Top: 
coupled (c) system. Middle: isolated (i) system where the 
couplings between the middle region and the leads are set to 
zero. Bottom: decoupled (d) system where the decoupling is 
achieved by introducing the potentials V£ and Vr. 



II. METHOD 

A. Non-equilibrium spin density 

The non-equilibrium Green's function (NEGF) ap- 
proach is based on dividing the junction into three re- 
gions (see Fig. [T]): two semi-infinite unperturbed leads 
left (L) and right (R) and a middle (M) region (or scat- 
tering region). This division allows a description of the 
effect on the middle region of the semi-inifinite leads each 
having a different chemical potential. The effect on the 
middle region can be written in terms of a self-energy of 
the left lead given by the coupling from the middle 
to the left lead and back 
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where gi, is the surface Green's function of the isolated 
semi- infinite left lead and vml describes the coupling of 
the left lead to the middle region. In an analogous way 
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one defines the self-energy of the right lead E#. These 
self-energies can be interpreted as fluxes of incoming and 
outgoing electrons at the connection between leads and 
middle region*! Using the self-energy of the left and right 
lead one can express the spin density matrix in the middle 
region 
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where Gm,c is the Green's function of region M cou- 
pled (c) to the semi-infinite leads and Jl and fn are 
the the distribution functions of the left and right lead»£ 
All quantities in Eq. [5] are energy dependent. The self- 
energy can be used to relate the coupled (c) and isolated 
(i) Green's functions of the middle region via a Dyson 
equation 
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This equation is based on the assumption that the self- 
energies can be added, which is true if the leads are well 
separated and their interaction can be neglected. 

Our approach for calculating the self-energy with 
a Korringa Kohn Rostoker (KKR) multiple scattering 
scheme is based on the approach due to Henk et al&, 
following the work by Pcndry et al£> Here we sketch this 
idea and highlight the basic assumptions which are nec- 
essary. 

For the calculation of the self-energy using Eq. © , one 
needs to know the coupling between the lead and the mid- 
dle region. In contrast to, e.g., tight binding approaches^ 
this coupling is not directly accessible within the KKR 
scheme because one calculates the Green's function of 
a system by a Dyson equation. Hence to obtain the 
coupling, one has to invert the Green's function. The 
alternate approach we use is to introduce decoupling po- 
tentials V d and V£ (see Fig. [T] bottom) which decouple 
the middle region from the leads using finite barriers. In 
the following we show that one can define a self-energy 
using V d and the Green's function Gd of the decoupled 
system. The Green's functions of the infinite systems 
(coupled and decoupled) can be related by the Dyson 
equation 



G c = G d - G d (V T d + V£) G c 



(4) 



Inserting this equation once in itself and using the as- 
sumptions (written schematically) 
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and 



G d (VtGdVt + V£G d V%) G c » 
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one can identify the self-energies by comparing the result 
to Eqs. © and © 



£i = y£G d Vt and S fl = v£ G d V&. 
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Assumptions © and © are necessary because V d and 
Vft are local potentials whereas vlm in Eq. © is a cou- 



pling. Assumption © is fulfilled if (V/ +V$) G c » 1 
which one assures by choosing an appropriately high V d . 
Assumption © is that the self-energies of the right and 
left lead can be added. This is fulfilled if the leads are 
well separated because the elements of G d relating the 
left and the right lead decay exponentially with respect 
to the thickness of the decoupling potential. 

By comparing Eqs. © and © the role of G d is the 
role of the surface Green's function of the isolated leads. 
Therefore, one can also set the whole middle region to 
the potential V d when calculating G d . For details of 
the implementation using the KKR basis set see Ref. H|. 
For the present work, we generalize the method to non- 
collinear magnetizations based on Ref. 1 1 Ot 



B. Spin-transfer torque 

From the non-equilibrium Green's function for a non- 
collinear magnetization described in the previous section, 
it is straightforward to compute the spin-tranfer torque. 
In linear response, the spin torque f per current / on 
layer i can be expressed by (see Haney et al£-) 
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where A; are the matrix elements of the exchange field 
along the magnetization axis of the layer i expanded 
into spherical harmonics with I being the angular mo- 
mentum. to' 1 * is the magnetic moment of the electrons 
contributing to the transport and is calculated from the 
non-equilibrium spin density matrix at the Fermi level 
E F using Eq. © 

P tr = 7^G C (E F ) (X L (E F ) - G\{E F ) (9) 

taking into account that the electronic states at Ep are 
occupied only in one lead. The transmission probability 
T is calculated by 



T = Tr 
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where the trace is over the spin index and the basis set 
expansion. Due to the in-plane translational invariancc 
of the junctions one can label the states by the wave 
vector ky and all quantities in Eqs. © and (TTO)) depend 
on k|| . To get the total values one has to integrate over 
the 2D Brillouin zone. 



III. APPLICATION TO Co/Cu/Co 

In this section we test our approach by applying it to a 
Co/Cu/Co system. In particular, we consider the depen- 
dence of the torque on the angle between the two mag- 
netizations of the ferromagnetic leads. For this purpose, 
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FIG. 2: Conductance of the Co/Cu/Co/Cu spin valve as a 
function of the relative angle between both magnetization of 
the ferromagnetic layers. 
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FIG. 3: Torque per current as a function of the relative angle 
between both magnetizations of the ferromagnetic layers for 
two different numbers of k)| points in the whole 2D Brillouin 
zone. 



we use the same structure used in Ref. [H consisting of 
a semi-infinite Co, 9 monolayers Cu, 15 monolayers Co, 
and semi-infinite Cu. The lattice constant is 0.361 nm. 
Fig. [U shows the dependence of the conductance g on the 
relative angle 9 between the magnetizations of the Co 
layers. We find excellent agreement with Ref. H and a 
GMR ratio (g(0°) - g(180°)) /g(180°) = 47 %. 

Fig.[3]shows the in-plane torque tii and the out-of-plane 
torque t_j_ as a function of the angle 9 for two different k|| 
point samples. The out-of-plane torque has two contribu- 
tions: one from the right going states that are occupied 
and have no left going counterparts, and the other from 
states below both chemical potentials in which both left 
and right going states are occupied. The latter contribu- 
tion requires integration over energy as well as parallel 
wave vector. However, the usage of a complex energy 
contour makes it easier to converge. For the former con- 
tribution one can not use a complex energy because only 
right going states are occupied. Therefore, to test our 
method against previous calculations and to test the k- 
point convergence, we consider only this contribution to 
the out-of-plane torque in the following. The dependence 
of tii on 9 is the same for both k|| samplings and is in 
good agreement to Ref. [H. In contrast, t± shows rapid 
oscillations as a function of 9 for the kii sampling using 
40 000 k|| points. A similar dependence was observed in 
Ref. 0. However, a significantly larger ky point sample 
leads to an almost smooth dependence of t± on 9. 

Due to the drastic change in tj_ from increasing the 
number of ky points, we test the convergence of the in- 
plane torque tii and the out-of-plane torque t± for a fixed 
angle 9 = 60° as a function of the number of k|| points 
(see Fig. For tm the convergence is fast and a rela- 
tively low number of k|| points is sufficient. On the other 
hand t± is very sensitive to the number of k|| points and 
a large number is necessary to get convergence. The slow 
convergence results from the presence of short period os- 
cillations at the corners of the 2D Brillouin zone which 
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FIG. 4: Convergence test for the torque with respect to the 
number of ky points in the whole 2D Brillouin zone. 



require a very fine mesh. Therefore, the rapid oscillations 
in t± as a function of angle found in Ref. [5| disappear for 
fully converged ky point samples. 



IV. CONCLUSION 

We present a method to calculate the spin-transfer 
torque within a screened KKR scheme by calculating 
the non-equilibrium spin density using the steady state 
Keldysh approach. The in-plane torque in the Co/Cu/Co 
junctions is robust with respect to the k|| point sampling 
but the out-of-plane torque converges slowly with respect 
to the number of ky points. The reason is that there are 
short period oscillations at the edges of the 2D Brillouin 
zone. These oscillation require a very fine ky point mesh 
to get the correct value for the integral. Then both com- 
ponents of the torque are a smooth function of the angle 
but with maxima at different angles. 
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